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In framework of statistical models, different particle ratios at energies ranging from 3.5 to 
200 GeV are calculated. Assuming that the particle production takes place along the freeze- 
out curve, we study the sharp peak in K + /n + ratio observed at SPS energy. We study the 
responsibility of non-equilibrium quark occupancy of phase space 7; for particle production. 
Allowing 7; to take values other than that of equilibrium, we got a very well description for 
a : K + /tv + ratio at all energies. Using the resulting parameter set, we analyzed the K /ty , 

A/tv + and A/ < ty > ratios. We found that the corresponding peaks all are located at the 
same value of energy, ^snn — 7.5 GeV. At this energy, the entropy per particle is singular. 
The saddle-point in entropy per particle likely refers to critical phenomenon and change in 
the phase space. 

m ; 

§1. Introduction 

^1- . 

'xf ■ 

£SJ ■ The statistical models have been successfully used to describe the particle pro- 

duction in different heavy-ion collisions. Studying the ratios of particle yields 1 )' 2 ) 

^ \ is of great interest, not only to determine the freeze-out parameters, and /xb, 

but also to eliminate the volume fluctuations and the dependence of the freeze-out 
surface on the initial conditions. On the other hand, serious challenges rise when 
Q_i| trying to bring together results from different experiments, like the dependence of 

T c h on fiB- Although each set of T^-^ib originally has been calculated by statistical 
models through combining various particle ratios comparable to the experimental 

^ \ results, it is still debated on the models which try to describe T c h vs. fis-^ 8 ^ 

The K + /ir + ratio represents another challenge for statistical models. The pre- 
diction of a sharp peak in K + /ir + ratio at SPS energy 9 )' 10 ) wasn't compatible with 
the statistical models. 11 )' 12 ) The latter merely expect a mild maximum, which ends 
up with a plateau. There were many attempts to interpret this phenomenon. 13 )' 14 ) 
In this work, we study the ratios of strangeness to non-strangeness particles. 
We restrict the discussion to the kaon-to-pion and lambda-to-pion ratios. In the 
experimentally estimated date, we observe that meanwhile there is a sharp peak 
in K + /n + ratio, the peak of K~/ir~ is much smooth. This is an indication of 
strangeness asymmetry, which should reflect itself in strange hyperon production. 
Therefore, a sharp peak in A/n + ratio is expected. The observed peak in Aj < tt > 
is smooth. Here < it >= 1.5(7r + + 7r~). References of experimental data can be 
found in 14 ) 

The experimental results on K + /tt + ratio at different energies will be fitted by 
7i. ji are the quark phase space occupancy parameters. They are allowed to take 
values other than unity; the equilibrium value. The subscript % refers to the light 
and strange quark flavors. This idea has been discussed by Johann Rafelski and his 
collaborators and implemented in their statistical hadronization model. 13 )' 15 ) Note 
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the differences between present work and Ref. 13 ^' 15 ^! We explicitly assume that the 
particle production takes place along the freeze-out curve, which is characterized 
by s/T 3 = 7, e ^ 8 ^ where s is the entropy density. This condition assumes constant 
degrees of freedom along the line of freeze-out, Fig.^ Using the parameter set from 
fitting K + /ir + vs. y/s, we calculate the other ratios. We find that the peak in 
K + /ir + is not a unique phenomenon. 9 )' 10 ) Its height and sharpness are obviously 
greater than that of the peaks observed in K~ /tt~, A/ir + and Aj < ir >. Another 
worthwhile finding is that all peaks are located at almost one value of energy. This 
might be an indication for a critical phenomenon. 

In the contrast to statistical models with ji = 1, our calculation with variable 
ji results in an excellent agreement with the particle ratios. For non-equilibrium 
freeze-out, i.e. variable j{, we find that the entropy per particle s/n which measures 
the averaged phase space density 16 ' 20 ' has a maximum value located at the same 
collision energy as the peaks of particle ratios do. This can be interpreted as a 
manifestation of critical phenomenon. According to recent lattice QCD calculations, 
the critical endpoint might be located at \ib ~ 0.4 GeV. The corresponding yfs 
is much close to the energy of the particle ratio peaks. Our estimation for s/n at 
RHIC energy is qualitatively consistent with the results given in 

Ref .18),21) Xo 

our 

knowledge, there is no experimental estimation for s/n at lower energies. But there 
is an indication reported in Ref. 18 ^ 21 - 1 that s/n at SPS is higher than at RHIC. This 
agrees very well with our predictions, Fig. [5J 
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Fig. 1. The freeze-out parameters T c h and calculated according to the condition of constant 
s/T s . In these calculations, 7, = 7^ = 1.0. 

§2. The model 

The pressure in hadronic phase is given by the contributions of all hadron res- 
onances treated as a free gas. 22 ) -25 -* At finite temperature T, strangeness and 
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iso-spin and baryo-chemical potential the pressure of one hadron reads 

9 f°° 2 T -e(fc) "I 

p(T, fjL B , US, Vl 3 ) = J k dkln\l± \ B \sXi s e T J, (2-1) 

where e(k) = (k 2 + m 2 ) 1 / 2 is the single-particle energy and ± stands for bosons 
and fermions, respectively, g is the spin-isospin degeneracy factor. 7 are the quark 
phase space occupancy parameters. A = exp(/z/T) is the fugacity. fj, is the chemical 
potential multiplied by the corresponding change. The particle number is given by 
the derivative of Ea. l2-ll with respect to the chemical potential of interest. The total 
pressure is obtained by summation all hadron resonances. 

The quark chemistry is given by relating the hadronic chemical potentials and 
7 to the quark constituents. 7 = 7™ 7™ with n and m being the number of light and 
strange quarks, respectively. fiB = 3// g and us = ^ q — fJ-s, with q and s are the light 
and strange quark quantum number, respectively. The baryo-chemical potential for 
the light quarks is fi q = (fi u + /id)/2. fis is calculated as a function of T and \xb 
under the condition of strangeness conservation. The iso-spin chemical potential 

A*/a = (Pu ~ Vd)l 2 - 

In carrying out our calculations, we use full grand-canonical statistical set of the 
thermodynamic parameters. Corrections due to van der Waals repulsive interactions 
has not been taken into account. 25 - 1 Although we do not use Boltzmann approxi- 
mation, we can for simplicity assume it, in order show on which parameters are the 
particle ratios depending. For finite iso-spin fugacity A/ 3 we get 

n K + _ K+ ! X q 



n n + 7T+ s Xj 3 
n-K- = 

n A A 



kK't*- - (2-2) 
ocA s f^) 21 (2 . 3 ) 

<* M -^-1 % 2 7 S (2-4) 



The particle numbers at zero chemical potential, rij(T) ~ Tm'jK2(mj/T) 1 represent 
the proportional factors in these expressions. nj(T) is a smooth function of T. The 
fugacity A is also smooth function of T. Correspondingly a monotonic increase in 
the particle ratios is expected with energy. 

As given in Ea. l2-H the statistical parameter 7 appears in the front of Boltzmann 
exponential, exp(— e/T). It gives the averaged occupancy of the phase space relative 
to equilibrium limit. Therefore, in the equilibrium limit 7 = 1. Assuming the time 
evolution of the system, we can describe 7$ as the ratio between the change in particle 
number before and after the chemical freeze-out, i.e. ji = rii(t) j rej(oo) . The chemical 
freeze-out is defined as a time scale, at which there is no longer particle production 
and the collisions become entirely elastic. In case of phase transition, 7^ is expected 
to be larger than one, because of the large degrees of freedom, weak coupling and 
expanding phase space above the critical temperature to quark-gluon plasma. 

As discussed above, thermal models using 7$ = 1 can't describe the characterized 
peak in K + /tt + . We display this in Fig. [3 
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Fig. 2. The K+/ir+, K~ /w~ , A/w + and A/ < tt > ratios at AGS (^5^7 < 4.84 GeV), SPS (6.26 < 
< 17.27 GeV) and RHIC energies (62.4 < ■ s /snn < 200 GeV). Both light and strange 
quark occupancy parameters 7 9 and 7 3 are assigned to unity. 

Allowing 7j to take values other than the unity makes it possible to relate ji 
to the collision energy. Then we calculate K + /tt + , Eq. (|2-2j) . The obtained values 
of 7g and "y s partly given in last two columns of Tab. [IJ are greater than one. This 
is an indication of oversaturation. To be concrete, we emphasize that 7q/7 s ratio 
first increases up to the collision energy at which the sharp peak of K + /tt + ratio is 
observed. Then it sharply decreases. For higher energies, it smoothly increases. The 
latter indicates to the fact that the mechanism responsible for the particle production 
remains unchanged from SPS to RHIC energies. In calculating the other particle 
ratios, we use the same parameter set. For K~ /n~ ratio, the ratio j s /7q multiplied by 
X s /X q , Eq. (|2-3j) . is the leading parameter. For A/ir + ratio, the occupancy parameters 
7g7 s , Eq. (|2-4j) . are responsible for the good agreement. 
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3.5 


0.8 ±0.05 


0.25 ±0.05 


0.9 ±0.02 


7.5 


1.0 ±0.10 


0.48 =F 0.07 


1.32 ±0.12 


17 


0.7 ±0.07 


1.6 =F 0.10 


1.43 ±0.13 


130 


0.8 ±0.08 


1.6 ±0.12 


1.48 ±0.11 



Table I. The values of 7 9 and 7 S at certain collision energies. The second column is related to the 
results depicted in Fig. [2] The corresponding 7 9 is constant and equals to one. The last two 
columns are the fit parameters related to Fig. 0] 

§3. Results 

In Fig. [2J we plot the predictions from statistical models, i.e. j q = j s = 1, on the 
top of the experimentally estimated results at AGS (y^ < 4.84 GeV), SPS (6.26 < 
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y/s < 17.27 GeV) and RHIC (62.4 < < 200 GeV) energies. For the K+/ir+ ratio, 
there is a good agreement at AGS energy. At the SPS energy, y/s ~ 10 GeV, there is 
a very mild maximum., which sets on to a plateau. At higher energies, this version of 
statistical models obviously overestimates the K + /ir + ratio*). 14 )' 27 ) Even by taking 
into consideration the error bars, we find that the resulting values are still above 
the "data". The same behavior can be seen in K~ /tt~ ratio (20% overestimation at 
y / sjv r /v = 130 GeV but still within the error bars!). The overestimation here starts 
earlier than K + /ir + . For strange baryons, we find a decreasing population of strange 
quarks with increasing energy. For A/ir + , there are underpredictions in two energy 
regions. The first one is at SPS energy. This indicates to a peak corresponding 
to that of K + /tt + . The second one is at RHIC energy (25% underestimation at 
y / s/v r /v = 130 GeV, for instance). Here the resulting values lie beneath the error 
bars. The same behavior is to be observed in Aj < ir >. 

We can so far summarize this discussion. The maxima do not all occur at the 
same energy and therefore, the case for a phase transition is not very strong. 1 ^ Fur- 
thermore, we conclude that the particle ratios calculated by statistical models with 
7g = 7s = lj disagree with the experimentally estimated results, especially at RHIC. 

In Fig. |31 we depict the results for variable "y s and constant 7 9 = 1. This as- 
sumes that only strangeness is out of equilibrium. The light quark numbers are in 
equilibrium. Using variable 7 S and j„ = 1 has a long tradition. 13 )' 28 )' 29 ) This has 
been based on assuming strangeness saturation as a signature for QGP. 30 ) Here, we 
want to establish the same procedure. We use ■js as a fit-parameter for K + /ir + 
ratio. The results are partly given in Tab. [H We find that 7 S always smaller than 
one. This means that the strangeness is always undersaturated. We will leave open 
the discussion on the physical interpretation, why the strangeness quantum number 
has to be out of equilibrium, while the light quark quantum numbers are explicitly 
in equilibrium? The results are convincing. An obvious underestimation for the 
particle ratios A/ir + and Aj < ir > is obtained. 

In Fig. 01 we draw the results with varying j{. In estimating the 7j-values, we 
fit the experimental results on K + /tt + ratio. j q is given as an input depending on 
T 13 )' 31 ) and correspondingly on tJsnn- We assumed that the particle production 
takes place along the freeze-out line, Fig. ^ at which we calculate the temperature, 
Tch(l^B), according to the condition that s/T 3 = 7. 6 )~ 8 ) s is the entropy density. All 
thermodynamic parameters at fixed. 7, are varying. 

It is interesting to note the overall agreement between our predictions and all the 
experimentally estimated particle ratios. The parameter set obtained from fitting 
K + /V + ratio are used to calculate the other particle ratios. There is no fitting done 
here. Although K~ /tt~ ratio raises with the energy, we notice a weak slope within 
the energy region 7 < ^/snn < 12 GeV. This might be an indication for strangeness 



*' For instance, at ^/sjvjv = 130 GeV, there is a 20% overestimation. This is in contradiction 
with the expectations given in Ref. 26) The value of K + /tt + ratio can be indirectly calculated from 
the results given in. 26 ' Accordingly, we find that K + /tt + ~ 0.f63 at the same energy. 
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Fig. 3. K + /ir + , K~ /n~ , A/n + and A/ < n > particle ratios as function of ^/sjvjv- 7 9 is assigned 
to one, meanwhile 7 S is a free parameter. Although the K + /ir + ratio has been successfully 
fitted, the calculation of the other particle ratios underestimates the experimental results. 

asymmetry. The same behavior can be seen in the experimental data as well. While, 
our calculations on this particle ratio still lie above the SPS results, the agreement 
with the AGS and RHIC results is excellent. We also notice that the K~/tt~ ratio 
is not so much sensitive to the change in the light quark occupancy parameter from 
■jq = 1 in Fig. |21 to the non-equilibrium values in this figure. 

The worthwhile finding here is that although the heights of the peaks are differ- 
ent (different strangeness asymmetries), all peaks are located at almost one value of 
energy, ^s^ (c) ~ 7.5 GeV. This energy value is corresponding to baryo-chemical 
potential of hb — 0.42 GeV. 14 ) This value has been reached by the lead beam 
at 40 AGeV at CERN-SPS. Another conclusion we can make so far is that non- 
equilibrium processes in both light and strange quark numbers are responsible for 
particle yields. 32 ) 

In present work, we study the responsibility of non-equilibrium quark occupancy 
of phase space for particle production. We devote this work to strangeness/non- 
strangeness particle ratios. In the following, we go one step farther in clarifying the 
physical reason behind the location of the peaks. The dependence of single-particle 
entropy on the collision energy is related to the averaged phase space density. In 
Boltzmann limit and for one particle 



where e is the single-particle energy. In this expression, s/n is related to e. Ap- 
parently, s/n becomes maximum when the chemical potential gets as large as the 
single-particle energy e. As we assumed Boltzmann limit, the maximum is unity in 
this case. Depending on the chemical potential fx, we can insert particles into the 
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Fig. 4. The K + /ir + , K~ /ir~ , A/ir + and A/ < n > ratios at energies ranging from 3.5 to 200 GeV. 
The experimentally estimated K + /tt + ratios have been fitted by y q and j s . The other particle 
ratios have been calculated by applying the resulting parameter set. 

phase space. The maximum occupation is reached at /i = e. This is an upper limit. 
Then it becomes prohibited to insert more particles. On the other hand, we can 
expect - at least theoretically - occupation values larger than this classical upper 
limit, only if the phase space itself is changed. This situation is most likely provoked 
by the phase transition. From this discussion we can apparently realize that s/n 
might play the same role as the statistical parameter 7 does. 

The results on s/n vs. y / ijv r /v are depicted in Fig. Again we use here many 
hadron resonances and full grand-canonical statistical set of the thermodynamic pa- 
rameters. In this case, the complete dependence of s/n on T and [i and consequently 
on y^s/vTv, can straightforwardly be obtained by deriving s and n from Eq. 12-11 For 
7 9 = 7 S = 1, we find that s/n increases as the energy raises from AGS to low SPS 
(\/ s nn < 10 GeV). As we saw in Fig. |2J a mild maximum in K + /ir + is located 
at the same value of ^/snn- For higher energies, s/n remains constant. Again, 
this also is the case in Fig. |2J This behavior might be an indication that a strong 
compensation of the collision energy takes place in this energy region. Although we 
introduce more energy in the system, the number of particles allowed to occupy the 
phase space remains constant. This is an indirect signature that the phase space 
itself remains constant. 

For the non-equilibrium case, i.e. variable j q and 7 S , we find almost the same 
behavior up to ^snn ~ 6.5 GeV. At this energy, there is a singularity in s/n. The 
singularity might be related to a certain critical phenomenon. 33 ) For the rest of SPS 
energies, s/n rapidly decreases. Although the energy is increased and consequently 
the produced particles, the single-particle entropy decreases. This can only be un- 
derstood, by assuming that the phase space shrinks. At RHIC, s/n decreases slowly 
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with the energy. The shrinking in phase space becomes slow. If this model would 
give the correct description, we now might have for the first time a theoretical expla- 
nation for the dependence of phase space on energy. The phase space at SPS energy 
is apparently larger than the phase space at RHIC and LHC. As discussed above, the 
same behaviour has been found, experimentally. 16 )~ 20 ) The consequences are that 
the quark-gluon plasma might be produced at SPS. And detecting its signatures at 
RHIC might be non-trivial. 

About the nature of critical phenomenon, we can't make any strong statement. 
One might think of phenomenological models 34 ) and lattice QCD calculations 35 )' 36 ) 
for critical endpoint. According to lattice, the endpoint so far might be at = 
0.36 ± 0.04 GeV. According to Ref., 36 ) \ib ~ 0.42 GeV. We have to mention here, 
that the lattice calculations exclusively assume equilibrium j q and j s . Therefore the 
comparison can't be straightforward. 
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Fig. 5. The entropy per particle s/n as function of ^/sjvjv- Only for variable quark occupancy 
parameters y q and 7 S there is a singularity in s/n ratio. It is important to notice that the 
singularity is located at almost the same energy as the peaks of particle ratios. 

§4. Summary and discussion 

We have studied the ratios of strangeness-to-pion in dependence on the collision 
energies from 3.5 to 200 GeV. We assumed that the particle production takes place 
along the freeze-out curve. Fixing all thermodynamic parameters and introducing 
quark occupancy parameters, we have fitted K + /ir + ratio. Using the resulting pa- 
rameters, we made predictions for other particle ratios. We found that the agreement 
is excellent. The S/ir and Q/ir ratios are shown in Ref. 13 ) Within the statistical 
acceptance, we conclude that almost all peaks are located at one value of energy. 

Different models have been suggested in order to interpret the maximum ob- 
served in the particle ratios. Besides the statistical model at early stage 9 )' 10 ) there 



Particle Ratios in Heavy-Ion Collisions 



9 



are two additional models. The first one relates the peak to a transition from baryon- 
rich to meson-rich hadronic matter. 14 ) According to the second model, 13 ) the peak 
separates a high entropy phase from a low entropy one. 

As shown, 7 g and 7 S played the ace in reproducing the particle ratios at all ener- 
gies. This implies that the particle production is due to non-equilibrium processes in 
both light and strange quarks. Assuming equilibrium in only one of these quantum 
numbers is not able to produce the particle ratios as given in Fig. 01 

We introduced s/n as a thermodynamic quantity related to the statistical pa- 
rameters 7. In Fig. [51 we found that the sharp peaks of the particle ratios are 
associated with a singularity in the entropy per particle. The latter can be under- 
stood by assuming a rapid modification in the phase space. The entropy per particle 
seems to play the same role as the statistical parameter 7 does. The dependence of 
phase space on energy is essential to access the phase transition. 
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